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Abstract 

The very recently observed D*j(2463) meson is described as a J p = 1 + cs bound state in 
a unitarised meson model, owing its existence to the strong OZI-allowed 3 Pq coupling to the 
nearby S-wave D*K threshold. By the same non-perturbative mechanism, the narrow axial- 
vector Z) s i(2536) resonance shows up as a quasi-bound-state partner embedded in the D*K 
continuum. With the same model and parameters, it is also shown that the preliminary 
broad 1 + Z?i(2400) resonance and the established narrow 1 + Di(2420) may be similar cn 
partners, as a result of the strong OZI-allowed 3 Pq coupling to the nearby S'-wave D*ir 
threshold. The continuum bound states Di(2420) and _D s i(2536) are found to be mixtures 
of 33% 3 Pi and 67% ^Pi, whereas their partners L>i(2400) and D* sJ (2463) have more or 
less the opposite +1 Pi-state content, but additionally with some D*ir or D*K admixture, 
respectively. 

The employed mechanism also reproduces the ratio of the K^-Ks mass difference and 
the Ks width, by describing Kl as a bound state embedded in the -kit continuum. The 
model's results for J p = 1 + states containing one b quark are also discussed. 

1 Introduction 

Bound states embedded in the continuum have been suggested by von Neumann and Wigner [1]. 
Since then many works on this theme have appeared in various fields of physics [2-6]. In the 
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present paper, we shall study such states appearing as (approximate) solutions of our simple 
unitarised meson model [7]. As three concrete applications, we choose here the K L -K S system, the 
-D s i(2536) [8] together with the very recently observed narrow L>*j(2463) [9,10] state, and finally 
the couple consisting of the established -Di(2420) [8] and the preliminary broad Di(2A00) [11,12] 
resonance. At first sight, it may seem odd to try to relate such utterly disparate states, ranging 
from mesons that can only decay weakly to very broad mesonic resonances. Moreover, while 
in the J p = 1 + cn (n stands for non-strange) system the _Di(2400) is much broader than the 
slightly heavier D 1 (2420), for the cs states the likely 1 + _D*j(2463) is even narrower than the 1 + 
/} s i(2536). Nevertheless, we shall demonstrate below that, in all three simple mechanism 

of coupling two degenerate qq channels with the same quantum numbers to the meson-meson 
continuum is capable of accounting for the experimental data, through an exact or approximate 
decoupling of one of the physical qq states. 

2 Degeneracy lifting via coupling to the continuum 

When two degenerate discrete channels are coupled to the same continuum channel, the de- 
generacy is lifted. One state decouples, partly or completely, depending on the details of the 
interaction dynamics, and turns into a continuum (approximate) bound state. The other state 
turns into either a resonance structure in the continuum, or a bound state below threshold when 
threshold is near enough. Within our unitarisation scheme, quark-antiquark channels are coupled 
to the meson-meson continuum by OZI-allowed 3 Pq qq pair creation and annihilation. 

In the present investigation, we confine our attention to two degenerate quark-antiquark 
systems, which can be distinguished by some internal structure irrelevant for the coupling to the 
meson-meson continuum. For example, the constituent quark masses of u and d are the same in 
most models. Hence, the internal dynamics of uu and dd vector states can in meson models be 
described by the same Hamiltonian. Under J PC = ++ quark-pair creation, both systems couple 
with the same intensity to pion pairs in P-wave. Nevertheless, because of the relative sign of the 
coupling constants under particle interchange [13], the (uu + dd) /y/2 state decouples from the 
two-pion continuum (uj meson), whereas the (uu — dd)/ y/2 state decays strongly into pion pairs 
(p° meson). 

Let us denote by ipi and ip2 the wave functions of the quark-antiquark systems, and by 
Hi and if 2 the Hamiltonians describing their internal confinement dynamics (e.g. employing a 
confinement potential). For the continuum we write the wave function ip and the dynamics H. 
The coupling interactions of the two confinement channels to the continuum are denoted by Vi 
and V 2 . In Ref. [13] it is found how the full spin, isospin and color degrees of freedom contribute 
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to the determination of such potentials. In the spirit of our model, we then obtain for the three 
channels the following set of coupled dynamical equations [14-17]. 

(H - E) i) + Vi + V 2 i)2 = 

V} i> + (#i - E) fa = 

V% i> + (H 2 - E) fa = . (1) 

We assume that the internal dynamics of the two degenerate quark- ant iquark systems does not 
depend on the difference in their internal structure, and, moreover, that these systems also couple 
the same way to the continuum, i.e., 

H 1 = H 2 and V x = a V , V 2 = (3 V (a 2 + (3 2 = 1) . (2) 

In this case it is easy to sub-diagonalise the system (1), so as to obtain 

(H-E)i> + V (c^i + PiJ 2 ) = 

(3) 

V* V + {Hi - E) (c^i + Pfa) = 

and 

(ifi - E) (^i - c# 2 ) = . (4) 

We end up with a system of the continuum ip coupled to a linear combination aipi + (3ip 2 of 
the confinement states in Eq. (3), and with a completely decoupled system for the orthogonal 
linear combination f3ip\ — aip 2 in Eq. (4). Since, moreover, Hi and H 2 are supposed to de- 
scribe confinement, Eq. (4) has only bound-state solutions, all embedded in the meson-meson 
continuum. 

3 The neutral kaon system 

A typical example of the above- described phenomenon is the two-pion decay mode of the neutral 
kaon system. Both ds and sd couple weakly to nn via the process depicted in Fig. 1. All other 
decay modes of neutral Kaons can to lowest order be neglected, since they couple orders of 
magnitude weaker. 

Because of particle-antiparticle symmetry, we may assume a — f3 — l/y/2. The CP = — 1 
combination (ds — sctj /y/2 completely decouples and turns into a bound state embedded in the 
mi continuum. This combination thus represents the Kl meson, which, since CP violation 
is not contemplated in our model, has no coupling to the two-pion continuum. However, the 
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Figure 1: W exchange takes care of the flavour transitions triggering the strong process of quark- 
pair creation in weak K°/K° two-pion decay. 



CP = +1 combination yds+sdj /\/2 exists as a, actually extremely narrow, resonance in S-w&vc 
7T7T scattering, hence describing the K$ meson. 

We assume here that the decay process depicted in Fig. 1 is dominated by the strong OZI- 
allowed 3 Po quark-pair-creation mechanism, whereas W exchange merely functions as a trigger to 
this process, which basically only determines the decay probability. Hence, besides the smallness 
of the parity- violating K$ — > ttti coupling constant, the transition potential describes here the 
coupling of a uu state with J PC = ++ , i.e., cr-meson quantum numbers to two pions, similar to 
our unitarised description of scalar mesons studied in Refs. [14,18,19] (see also Ref. [20]). In 
the delta-shell approximation for the transition potential in the case of the scalar Kq(800) and 
a (980) resonances [18, 19], we obtained for the delta-shell radius a a value of about a = 3.2 
GeV -1 . We shall hold on to this value in the following. 

A general solution of Eq. (3) for the S'-wave scattering phase shift cotg(5(s)), as a function 
of the total invariant two-pion mass y/s, is in this approximation and for small coupling given by 



cotg (S(s)) 



E + R(s) 


^0 


2" 




I(s) 


l^ol 


2 



(5) 



where E Q represents the ground-state energy of the ds (sd) when uncoupled, hence the mass of 
the K L meson. The remaining factors R(s), I(s), and JF are well explained in Ref. [21]. 

From expression (5) it is easy to perturbatively extract the real and imaginary part of the 
resonance pole position in the complex-energy plane, i.e., 



£pole ^ E + AE , with $le(AE) = R(s) \F \ and ^rn(AE) = I(s) \F \ Z . (6) 

The position of the K s resonance pole with respect to the K L mass is shown in Fig. 2. Since the 
coupling of the neutral kaons to the two-pion continuum is extremely small, the arrow pointing 
from y/s = mi to the Ks resonance pole position represents the pole trajectory for increasing 
intensity of the transition potential Vi, too. 
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Figure 2: Kl bound-state pole embedded in the continuum, and Ks resonance pole in the second 
Riemann sheet. 

For the ratio of the Ks decay width Ts, equalling twice the imaginary part (6) of the resonance 
pole position in the complex-energy plane, and the neutral-kaon mass difference, which equals 
the real shift (6) of the pole with respect to E = m L , we read from formula (5) the result 



(7) 



m s - m L R(s) 

From Ref. [21] we moreover learn that, in the case of S'-wave scattering, one has 

M = 3 ^tA = -t§M • (8) 
R(s) n {ka) BK J K ' 

As a consequence of Eqs. (5), (7), and (8), we obtain an extremely simple relation for the mass 
difference in the neutral kaon system, the width of the Ks meson, the two-pion momentum k 
and the strong-interaction radius a, reading 



w tg(fca) . (9) 
m L - m s 

Substitution of k — 0.206 GeV and a = 3.2 GeV -1 gives us then the result 

-r 

— 2 — — w 1.3 (Experiment: 1.06 [8]) . (10) 
m L -m s 

Similar conclusions can be found in Refs. [22-26], where long-distance effects have been studied 
in more detail. 

4 The J p = l+ cs states 

In the case of the J p = 1 + cs states, we deal with two distinct systems, 3 Pi and 1 P 1 , which, as 
far as confinement is concerned, are degenerate when ignoring possible spin-orbit effects. Both 
states couple strongly to D*K, with threshold at about 2.501 GeV. 
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For the neutral kaon system, we could straightforwardly assume that the transition potentials 
V\ and V 2 in Eq. (1) are equal, because of particle-antiparticle symmetry. But for the 1 + cs 
systems we do not have such simple arguments. Nevertheless, if V\ and V2 are proportional, we 
still are in a situation comparable to the one discussed in the previous section, though now for 
strong interactions. 

In Ref. [13] it is shown that the spatial form of the transition potential mainly depends on 
the orbital angular momenta of the cs and D*K channels, whereas the relative intensities follow 
from recoupling coefficients. For transitions to vector+pseudoscalar we find that the intensity 
for 3 Pi is twice as large as for 1 P\. Consequently, by the use of Eqs. (4) and (3), we conclude 
that the J p = 1 + cs bound state embedded in the D*K continuum consists of a mixture of 33% 
3 Pi and 67% l P\, whereas its supposed resonance partner has more or less the opposite content, 
but additionally with some D*K admixture. 

In Sec. (3) we have used Eq. (3) for the description of tttt scattering in the presence of the 
weak coupling to the neutral kaon system. Here, we assume that Eq. (3) is also suited to describe 
D*K scattering in the presence of an infinity of cs confinement states. In our unitarised model, 
this just implies substituting the effective nonstrange quark mass by the charmed quark mass, 
the pion masses by the D* and K masses, and changing some of the quantum numbers. The 
scattering phase shift is given by an expression similar to the one shown in Eq. (5), but, since 
the interaction is not weak now, also higher radial excitations of the cs confinement spectrum 
must be included. Thus, we use the general expression [14, 18] 



cotg (5(s)) 



n (pa) 
Jo (pa) 



2\ 2 jipajl(pa) 



;n) 



n=0 E n 

Expression (11) is not only valid above threshold, where for small coupling (A <C 1) it gen- 
erates Breit-Wigner-like resonance structures in the scattering cross section around each of the 
energy eigenvalues E n (n = 0, 1, 2, . . .), but it is also valid below threshold. Resonances are 
characterized by complex-energy poles which are given by solutions of cotg(5(s)) = i. When 
solutions of cotg(<5(s)) = i come on the real-energy axis below threshold, then they describe the 
real energy eigenvalues for the bound state solutions of equation 3. For small values of A the 
latter poles are found in the proximity of those E n (n — 0, 1, 2, . . .) which come below threshold. 

The non-relativistic form of the Flatte formula [27] can be obtained from formula (11) in the 
case of overlapping resonances. But, then its analytic continuation to below threshold is lost. 

From expression (11) one can study the behavior of resonance poles, just being solutions of 
cotg(£(s)) = i, in function of variations of the model parameters (see e.g. Ref [21]). Their 
passage from above threshold, where resonance poles exist in the the second Riemann sheet, to 
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below threshold, where poles are expected to reside on the real-energy axis in the first Riemann 
sheet, is smooth and without any difficulties. For S'-wave scattering the ground-state pole hits 
the real-energy axis below threshold, then, for increasing values of A, moves along the real-energy 
axis towards threshold where it changes Riemann sheet, turns 180 degrees and starts moving 
towards smaller values of energy. 

Here, similarly to the procedure outlined in Ref. [14,18], we approximate the full sum over 
all cs confinement states by the two nearest states, that is, the ground state at E and the first 
radial excitation at E ± , plus a rest term. The latter is scaled to 1 by absorbing its value into the 
coupling constant A, yielding 



CO I *p CS I 
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Figure 3: Trajectory of the ground-state pole in the D*K scattering amplitude, as a function 
of the coupling constant A, defined in Eq. (11), and which parametrises the intensity of the 
transitions between the pure cs system and the D*K continuum (see Eq. 3). The arrow indicates 
how the pole moves when A is increased. For vanishing coupling, as well as for the confined 
system described by Eq. (4), the pole is on the real E axis. For large coupling we find the pole 
on the real E axis below threshold. Units are in GeV. 

We determine the ground-state mass of the uncoupled system (4) from the model parameters 
(m s = 0.508, m c = 1.562 and to = 0.19 GeV) given in Ref. [16], yielding E =2.545 GeV. This 
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also means that the bound state embedded in the D*K continuum has exactly the mass E , only 
9 MeV away from the experimental _D sl (2536) mass [8]. The first radial excitation lies 2lu = 380 
MeV higher. 

In Fig. 3 we study how the position of the ground-state singularity in the scattering amplitude 
for Eq. (3) varies, as the intensity of the transition potential V\ is changed. In Nature one 
can only "measure" the pole position for one particular value of this intensity, given by strong 
interactions. Nevertheless, it is very illustrative to study alternative values. In particular, we 
notice that for larger intensities of the transition mechanism the pole comes out on the real E 
axis, below the D*K threshold. Had we started from the perturbative formula (5), where R(s) 
and I(s) are proportional to A 2 , we would have started out as in Fig. 2, and never returned to 
the real E axis for increasing values of the transition intensity. The trajectory in Fig. 3 is highly 
nonperturbative, which can only be achieved when all orders [21] are accounted for. The masses 
obtained for the model's values of A are indicated by • in Fig. 3. For the _D*j (2463) we read from 
the figure m = 2.446 GeV, so, after all, this state is not a resonance as one would naively (i.e., 
perturbatively) expect from the coupled set of equations (3), in agreement with experiment. As 
mentioned above, for the .D s i(2536) we obtain 2.545 GeV from pure confinement. 

When the transition potentials V\ and Vi in Eq. (1) are not proportional, one has no simple 
diagonalisation, since commutators with the Hamiltonians will spoil the simplicity. Also, if Hi 
and H 2 in Eq. (1) are not equal, diagonalisation will not lead to completely decoupled systems. 
However, from experiment we learn that the D*K width of the -D s i(2536) is small (less than 
2.3 MeV [8]), implying that the bulk of the interaction indeed stems from the unitarisation 
mechanism. 



In order to describe the ground states of the J p = 1 + cn spectrum, we use the same Eqs. (5) 
and (6) as in the previous case. Also here, we determine the ground-state mass of the uncoupled 
system (4) from the model parameters (m U:d = 0.406, m c = 1.562 and uo = 0.19 GeV) in Ref. [16], 
which now yields E =2.443 GeV. Hence, the J p = 1 + cn bound state embedded in the D\ 
continuum comes out, in our model, at 2.443 GeV, some 20 MeV above the experimental [8,11] 
Di(2420) mass. The first radial excitation lies, as before, 2u = 380 MeV higher. The remaining 
parameters a and A are kept the same as in the previous case, yet scaled with the reduced 
constituent qq mass /i qq , in order to guarantee flavor invariance of strong interactions, i.e. 



5 The J p = 1+ cn states 




constant and 




constant 



(13) 
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where x and y represent the two flavors involved. The constants of formula (13) are fixed by [19] 
a us = 3.2 GeV- 1 and X us = 0.75 GeV" 3 / 2 for S wave. 

The resulting cross section is given in Fig. 4. Our peak shows up somewhat above 2.3 GeV, 
whereas the width of our D 1 (2420) is about 200 MeV. The corresponding experimental values are 
2400 ± 30 ± 20 MeV [11] (2461±i MeV [12]), and 380 ± 100 ± 100 MeV [11] (2901^° MeV [12]), 
respectively. In Fig. 4 we can further observe that the first radial excitation of the system 
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Figure 4: Model result for the cross section in units of GeV 2 for elastic S'-wave D*ir scattering, 
as a function of the total invariant mass E in units of GeV. 

of equations (3), for J p = 1 + cn, comes out more than 500 MeV higher than the ground-state 
resonance. Nevertheless, in expression (6) we always used Ei — E = 380 MeV [16]. Such a highly 
non-perturbative behaviour is inherent in the unitarisation procedure leading to formula (5). 

The bound state Z) 1 (2420) embedded in the D*ir continuum has zero width in our model, and 
so is invisible in the cross section of Fig. 4. On the other hand, experiment finds 18. MeV [8] 
and 26.7 ± 3.1 ± 2.2 [11] MeV for the full width, which is nonetheless very small as compared to 
the available phase space. Consequently, also in this case our assumption in Eq. (2) appears to 
be reasonable. 

The relative intensities a and (3 in Eq. (2) are the same as in the cs case. Hence, the relative 
content of the £>!(2420) is again 33% 3 Pi and 67% x Pi. The D^OO) roughly has the opposite 
mixture, plus a D*n component. 

In Fig. 5 we show how the scattering pole moves in the complex-energy plane when A is varied, 
which is the pole associated with the cross section of Fig. 4. 
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Figure 5: Trajectory of the ground-state pole in the D*tt scattering amplitude, as a function 
of the coupling constant A. The arrow indicates how the pole moves when A is increased. For 
vanishing coupling, as well as for the confined system described by Eq. (4), here representing 
the _Di(2420), the pole is on the real E axis. The pole postion for Eq. (3), which describes the 
.Di(2400), is indicated by • for the physical value of A. Units are in GeV. 

6 J p = l+ systems containing one b quark 

For J p = 1 + systems containing one b quark we may repeat the above described procedure. We 
determine the ground-state masses of the uncoupled systems (4) from the model parameters 
[m ud = 0.406, m s = 0.508, m c = 1.562, m b = 4.724 and u = 0.19 GeV) given in Ref. [16]. The 
remaining parameters a and A are kept the same as for the case of the D mesons. This results in 
5.707, 5.809 and 6.863 GeV for respectively bu/d, bs and be. The model's results for the partner 
states which couple to the continuum, are collected in table (1). Similar to the case of D*ir 
(section 5), we find in B*n a resonance structure above threshold, peak value at 5.6 GeV and 
some 120 MeV wide. The experimental result might be peaking at a slightly higher value [28-31], 
but the width of this resonance is in good agreement with our result. For the continuum bound 
state, which we find at 5.71 GeV, the OPAL collaboration reports a narrow resonance at 5.74 
GeV [31], whereas the CDF collaboration gives 5.719 GeV for the narrow B\ resonance [30]. We 
also might compare with the prediction of Kalashnikova and Nefediev, who obtain 5.716 and 
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5.741 GeV for the two B l states [32]. 

It should be noticed from table 1 that the 1 + bs partner state comes in our model at 5.73 
GeV, which is in the same energy region as the other two states. However, the CDF collaboration 
reports a very small (3.7%) contamination due to B s . 



system 


threshold 
GeV 


continuum 
bound state 

GeV 


partner state 
GeV 


scattering 
length 

GeV" 1 


bu 


5.465 (B*°n + ) 


5.71 


5.59 - 0.06i 


-1.13 


bs 


5.819 (B*-K+) 


5.81 


5.73 


2.08 


be 


7.190 (B*°D+) 


6.86 


6.83 


-31.8 



Table 1: Purely bound states and partners for J p = 1 + systems containing one b quark. 

Literally speaking, the 1 + B s and B c do not make part of the combinations continuum- 
-bound-state/partner-state, since even the purely bound states come below threshold. However, 
the case of bs could be different experimentally, given the error of some 10-30 MeV in our model 
predictions. So, there might be a continuum bound state, slightly above the B*K threshold. 
But, the partner state is definitely below threshold and thus very narrow. 

For B c both states are well below threshold. 

7 Summary and conclusions 

In the present paper, we have studied examples of mesonic systems that generate bound or quasi- 
bound states in the continuum, within an unitarised quark-meson model. In the neutral-kaon 
system, our approach provides a simple explanation for the widths and the mass difference of 
the K L and K$, which is also in agreement with the conclusions of more sophisticated methods. 
Application to the J p = l + cs and cn axial-vector charmed mesons accomplishes a simultaneous 
description of the established narrow _D s i(2536) and _Di(2420) states, as well as the recently 
observed very narrow D*j(2463), assuming it indeed is a 1 + state, and broad _Dx(2400), which is 
rather problematic in standard quark models. 

We predict the continuum bound states -Di(2420) and L> s i(2536) to be mixtures of 33% 3 Pi and 
67% ^i, which might be measured through radiative transitions [33]. For their partners _Di(2400) 
and .D*j(2463) we predict more or less the opposite 25+1 Pi-state content, but additionally with 
some D*n or D*K admixture, respectively. This is in agreement with the mixtures proposed by 
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Godfrey and Kokoski in Ref. [34], based on the flux-tube-breaking mechanism. Their 9 = —38° 
for cs amounts to 38% 3 Pi and 62% 1 P i regarding the state degenerate with 3 P 2 , and the opposite 
content for the state degenerate with 3 P . The former mixture corresponds to the continuum 
bound state P s i(2536), which is indeed approximately degenerate with the 3 P 2 in our model, since 
this J = 2 state is subject to only very small unitarisation effects due to the P-wave D*K channel. 
However, the latter mixture corresponds to the partner state D* sl at 2.443 GeV, which is strongly 
influenced by the coupling to the S'-wave D*K channel, just as the likely 3 P D*j(2317) + [35] 
is drastically affected by the S'-wave DK threshold [36]. Therefore, no approximate degeneracy 
holds for these two charmed mesons. Concerning the cu states, Godfrey and Kokoski obtained 
9 = —26°, resulting in 19% 3 Pi and 81% 1 P 1 , which indicates that our assumption (2) is probably 
somewhat less accurate in this case. 

For J p = 1 + quark- ant iquark systems which contain one b quark, we obtain good results as 
far as we can compare to experiment. 

In conclusion, we have once again demonstrated that unitarisation has to be incorporated 
in realistic quark models of mesons and baryons, as it constitutes the second most important 
interaction after confinement. 
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